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Dynamics of the 2d Potts model phase transition* 
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The dynamics of 2d Potts models, which are temperature driven through the phase transition using updating 
procedures in the Glauber universality class, is investigated. We present calculations of the hysteresis for the 
(internal) energy and for Fortuin-Kasteleyn clusters. The shape of the hysteresis is used to define finite volume 
estimators of physical observables, which can be used to study the approach to the infinite volume limit. We com- 
pare with equilibrium configurations and the preliminary indications are that the dynamics leads to considerable 
alterations of the statistical properties of the configurations studied. 



1. Introduction 

Lattice gauge theory investigations of the de- 
confining phase transition, for reviews see [Q, 
have mainly been limited to studies of their equi- 
librium properties. In contrast to that, their oc- 
currence in nature is governed by temperature 
driven dynamics. In the early universe it is a cool- 
ing process, which is estimated to proceed slowly 
on the scale set by hadronic relaxation times. In 
heavy ion collisions a rapid heating of the nuclei 
is followed by a much slower cooling process and 
it is no longer clear whether this cooling is slow 
or fast on the relaxation time scale. In addition, 
the small sizes of the nuclei systems impose com- 
plications which require to be studied. 

What dynamics should we study? Time evolu- 
tion in nature would correspond to some molecu- 
lar dynamics simulation of full QCD in Minkowski 
space, which is out of the question as already Eu- 
clidean space Monte Carlo (MC) calculations of 
QCD stretch computer resources to their limits. 
One may, however, argue that the computer-time 
evolution of a MC simulation reveals the correct 
physical features, as long as the updating proce- 
dure stays in the universality class of Glauber dy- 
namics [^. These are local MC updating schemes. 



which imitate the thermal fluctuations of nature. 
Examples are single- and multi-hit Metropolis al- 
gorithms, as well as the heat-bath algorithm. 

Also to simulate the dynamics of the decon- 
fining phase transition within Euclidean lattice 
QCD is out of reach for an exploratory study. To 
gain some qualitative understanding, we have de- 
cided to investigate g-state Potts models, whose 
simulations are far less computer time intensive. 
The relationship to QCD is that an effective 
model for the Polyakov loop is provided by the 
3(i, 3-state Potts model in an external magnetic 
field . Even this simplification is not yet a con- 
venient starting point. To get confidence in our 
numerical methods, we limit our first round of 
simulations to the g-state Potts models in 2c?, for 
which a number of rigorous results exist |^,|^ . 

In the following we drive our systems many 
times through various heating and cooling sched- 
ules and measure physical observables along the 
way. Averages and their error bars are calculated 
with respect to the repetitions of the schedule. 
Our observables are the internal energy and prop- 
erties of Fortuin-Kasteleyn |^ clusters. 

Details of our investigations are reported in the 
next section and a brief summary with conclu- 
sions is given in the final section. 



*This work was in part supported by the US Department 
of Energy under contract DE-FG02-97ER41022. 
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2. Results 

We have performed simulations on L x L lat- 
tices for the 2d Potts model with q = 2, 5 and 10, 
to have an example of a second order transition 
{q = 2), a weak {q = 5) and a strong {q — 10) 
first order transition. Due to space restrictions 
we limit the following presentation to g = 10. 

For suitably chosen values of /3min and /3max, we 
run in the range /3min < < /3max through each 
hysteresis cycle at least 640 times with f3 changed 
by ± A/3 after every sweep. The stepsize is 



A/3 



2(/3n 



/^min) 



nf3 



L2 



(1) 



where (3niin and Pniaic define the terminal temper- 
atures and the integer np = 1,2,... is varied. In 
the limit where we send L'^ —> oo first and then 
A/3 ^ the hysteresis is expected to survive, 
whereas it disappears if we send A/3 — > first 
and then — > oo. 
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Figure 1. Energy hysteresis curves for the 10- 
statc d ~ 2 Potts model on a 80 x 80 lattice for 
/3min = 0.4, /3niax = 1-2 (uot the entire range is 
shown) and the values of rt^ as given in the figure. 



For a first order phase transition, the slow- 
ing down of the Markov process is exponential 
in computer time, ~ exp[2/sL], where fs is the 
intcrfacial tension. In the following we study the 
limit L oo and A/3(L; np) — > for fixed. 
For 1-hit Metropolis updating we give in figure |l| 



our energy hysteresis data for the 10-state Potts 
model on a 80 X 80 lattice. Towards large np we 
see convergence to the equilibrium curve, which 
is obtained by means of a multicanonical sim- 
ulation. To analyze the physical content of such 
a hysteresis, we define for each np finite volume 
estimators for a physical observable as its value at 
the maximum opening of the corresponding hys- 
teresis. 
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Figure 2. Finite size latent heat estimates for the 
10-state d = 2 Potts model versus L for fixed val- 
ues of obtained from the maximum opening 
of energy hysteresis curves like those of figure |l|. 
The lines are just to guide the eyes. 



Figure ^ shows, for various fixed values of n^, 
the thus obtained latent heat estimates versus the 
lattice size L. For sufficiently large L all the esti- 
mators should, independently of their np values, 
converge towards the infinite volume value, which 
is also given in the figure. For each n/3 we have 
performed fits to extrapolate towards L ^ oo, 
but only some of them give good estimates. For 
our small and medium sized lattices the results 
appear to suffer from finite size corrections in L 
as well as in A/3, with the additional twist that 
A/3 can also be too small for some of the lattice 
sizes used. The best L —t oo estimate is 



AE = 0.3485 (12) at np = 16 . 



(2) 



Results for other values of q, other observables 
and details of these fits will be reported elsewhere. 
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To focus on a detailed physical understanding 
of the finite volume transitions, we have begun 
a cluster analysis of the configurations. The ap- 
propriate cluster definition is the statistical one 
of Fortuin and Kasteleyn j|] , see Ref. (D for a re- 
view of this and related topics. We are inter- 
ested in the effects of heating and cooling on the 
cluster structure, compared to their equilibrium 
structure. Whereas nucleation theory holds in the 
mean field approximation, phase transitions can 
also proceed via spinodal decomposition as the 
result of an off-equilibrium quench, which may 
still lead to observable signals, even when there 
is no longer a transition in the strict thermody- 
namical sense. This is somewhat different from 
the program of Satz ^ , who focuses on geomet- 
ric clusters and would like to use them as a signal 
when there is no proper transition. 




Figure 3. The probability of having a percolating 
cluster for the 10-state d — 2 Potts model on 
various lattice sizes (n^'' = 1). 



For the 10-state model we show in figure || our 
result for percolation of the ordered cluster. A 
heat bath algorithm was employed with n'^' — 1. 
The figure provides evidence that the geometry of 
the Fortuin-Kasteleyn clusters is distinct during 
the cooling and the heating part of the hysteresis. 
The percolation probability of the ordered phase 
still increases, when we already heat the system 
up again. A detailed investigation of the related 
cluster properties is in progress pO||. 



3. Summary and Conclusions 

Our energy hysteresis method allows for dy- 
namical estimates of the latent heat and works 
similarly for other physical observables. While 
the precision of these estimates is not competi- 
tive with those of equilibrium investigations 0, 
the hysteresis method provides information about 
deviations which are rooted in the dynamics of 
the heating and cooling process. From our ongo- 
ing analysis of the Fortuin-Kasteleyn clusters we 
gain insight into the dynamics of the transition, 
which may help to identify signals observable in 
heavy ion experiments. 
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